Long- and short-range correlations in the ab-initio no-core shell model 
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In the framework of the ab-imtio no-core shell model (NCSM), we describe the longitudinal- 
longitudinal distribution function, part of the inclusive (e, e') longitudinal response. In the two-body 
cluster approximation, we compute the effective operators consistent with the unitary transformation 
used to obtain the effective Hamiltonian. When short-range correlations are probed, the results 
display independence from the model space size and length scale. Long-range correlations are more 
difficult to model in the NCSM and they can be described only by increasing the model space or 
increasing the cluster size. In order to illustrate the model space independence for short-range 
observables, we present results for a large set of model spaces for *He, and in — AhQ. model spaces 
for 12 C. 

PACS numbers: 21.60.Cs, 23.20.-g, 23.20.Js 



Atomic nuclei are the result of a delicate interplay be- 
tween short- and long-range correlations among the nu- 
cleons, which makes their theoretical description rather 
challenging. For light nuclei, very successful methods, 
such as Green's function Monte Carlo 0, hyperspheri- 
cal harmonics 0], and the no-core shell model (NCSM) 
0, 0, 0, 0, 01 , have been developed recently. They allow 
an ab-initio description of nuclear properties, the only 
ingredients being realistic nucleon-nucleon (NN) interac- 
tions, which describe the experimental phaseshifts with 
high accuracy, and theoretical three- nuclcon forces. 



In the ab initio NCSM, one starts with a realistic 
NN interaction (theoretical three-body forces can also be 
used, but we will not discuss this case here) and per- 
forms a unitary transformation |^ 0| to a model 
space, which allows an exact diagonalization in a finite 
many-body space, defined by the number of excitations 
above a mean-field-like configuration. Details of the pro- 
cedure are available to the interested reader in previous 
publications 0, 0, IE Hi ■ Recently, we have extended the 
same procedure from the Hamiltonian to general one- 
and two-body operators 0, 0, 0| . The effect of the 
procedure is to reduce the dependence of the observables 
upon the model space and harmonic oscillator (HO) fre- 
quency, and, in the lowest approximation, it has proven 
to be effective only for short-range operators ^3 • 

Using the unitary transformation approach 0, 0, 0| , 
we obtain the following expression for the effective oper- 
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where the transformation operator uj satisfies the condi- 
tion QluP ~ to, with P and Q the projector operators in 
the model and complementary spaces, respectively, and 
O is the bare general operator, which acts in the entire 
space. 

In principle, the transformation operator uj can be 
computed using a finite set of eigenvectors in the full 
space, making use of the overlap of the full space eigen- 
vectors with the basis states in the P and Q spaces 0,0, 



I.e., 
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In the last equation, \ap) and {aq) are the basis states in 
the P and Q spaces, respectively. The summation runs 
over a finite subset, /C, of eigenvectors in the full space, 
and the tilde stands for the inverse of the overlap matrix, 

Equation (|2Jl shows that in order to obtain the trans- 
formation operator co one needs the solution to the initial 
A-body problem. This makes its application impractical, 
unless we use approximations. In the simplest approxi- 
mation, the transformation operator w and, therefore, 
the effective interaction are obtained in the relative sys- 
tem of two particles, in a large HO basis. The Q space 
is chosen to be a few hundred Ml, excitations in order 
to obtain an exact solution to the two-body Schrodingcr 
equation. Due to the rotational symmetry, we formu- 
late the problem in two-nucleon channels with good total 
spin s, total angular momentum j, and isospin t, reduc- 
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ing drastically the dimensions involved, when perform- 
ing the summation over the states in the Q space in Eq. 

The same procedure can be applied to operators 
that can be analytically expressed in terms of relative 
and center-of-mass coordinates of pairs. However, note 
that in the case of non-scalar operators calculations per- 
formed with Eq. become more difficult, because such 
operators can in general couple different channels. As ex- 
pected, since the transformation operator is a scalar, this 
procedure changes neither the character nor the rank of 
the bare operator. By construction, keeping the cluster 
approximation fixed (in this case, the two-body cluster) 
and increasing the model space decreases the effect of the 
renormalization. This was demonstrated in earlier papers 
in the case of ■^H for the ground-state energy 0, and in 
this paper we will also demonstrate this by comparisons 
between the effective and bare results for an observable 
related to the Coulomb sum rule. 

The inclusive (e, e') longitudinal data presents one of 
the clearest experimental signatures for short-range cor- 
relations in the wave-function of the ground state, at least 
for light nuclei. The Coulomb sum rule 



1 [°° 

Shiq) = — / duiSL{q,ui) 



(3) 



is the total integrated strength measured in electron scat- 
tering. InEq. lEl, S'l((7,w) = w)/|G£;,p((7, with 
R{q, Lo) the longitudinal response function and GE,p{q, lo) 
the proton electric form factor, while ujei is the energy 
of the recoiling A-nuclcon system with Z protons. The 
Coulomb sum rule SL{q), which is related to the Fourier 
transform of the proton-proton distribution function , 
can be expressed as [l6| 
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^(5-s-|pI(q)pL(q)|.g.s.) ^ 

= 1 + PLhiq) - ZFL{q)/GE,p{q,(^el) 
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where FE{q) is the longitudinal form factor. If one ne- 
glects the relativistic corrections and two-body currents, 
/9L(q) is simply the charge operator 



PL(q) 
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exp(iq • ri)(l -I- r^^i). 



Consequently, the longitudinal-longitudinal distribution 
function becomes 

PLhiq) = -7J^{9-S-\k{q\'Ci-Yj\){\^T^,{){l^T^^j)\g.s). 
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We present the results for PLL{q) for ^He in Figs. 
and El We have limited this investigation to two-body 
interactions only, and, in particular, have used the phe- 
nomenological CD-Bonn NN force because it yields 
reasonable convergence properties with increasing the 
size of the model space. Although experimental data for 
the longitudinal- longitudinal distribution function exist. 
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FIG. 1: (Color online) The longitudinal-longitudinal distri- 
bution function pLL(q) in *He for two model spaces (WhO, 
and 18hfl) and fixed frequency hfl — 28 MeV, using bare 
(dashed curves) and eflfective (continuous curves) operators. 
As discussed in the text, the results obtained with eflfective 
operators are almost indistinguishable. 



"1 ' \ ' I ' I : 

hQ = 19 MeV 

mi = 28 MeV; 



io;if2 
whn 
i8hn 
i8hn 




10 I 
10-' I 

10-^ 

I 10" 
10-' 
10-^ 
10-^ 
10-^ 
10-^^ 



FIG. 2: (Color online) Longitudinal-longitudinal distribution 
function, using bare (upper panel) and eflfective (lower panel) 
operators. We used two different HO frequencies, 19 MeV 
and 28 MeV, and two model spaces, lOhfl and 18hQ. 



a direct comparison with experiment is not suitable, be- 
cause we neglect: on one hand, (i) three-body forces in 
the model Hamiltonian, and, on the other hand, (ii) ex- 
change currents and (iii) relativistic corrections for the 
charge operator. Nevertheless, the results demonstrate 
the behavior of short- and long-range operators within 
the framework of the NCSM. 

Using a Gaussian operator of variable range, we have 
shown previously how the renormalization of this two- 
body operator depends upon its range We found 
that a short-range two-body operator is renormalized ac- 
curately at the two-body cluster level, while a long-range 
operator is weakly renormalized. The same behavior can 
also be inferred from Fig. ^ where we present PLhiq) 
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FIG. 3: (Color online) The convergence of the ground-state 
energy for the two frequencies used to compute the pll{<i)- 
The dotted line is the exact ground-state energy for the CD- 
Bonn interaction (-26.16 MeV p^l. 



calculated in two model spaces for a fixed frequency, us- 
ing both bare and effective operators. Thus, at large 
momentum transfer, the use of effective operators pro- 
duces model-space independent results. Even in small 
model spaces we obtain good results, although Fig. O 
shows that the ground-state wave function is not fully 
converged in such small spaces, since the ground-state 
enegy is not converged to the exact value. In particular, 
for Nraax = 10 (or lOftfi, in terms of allowed excitations 
beyond the lowest configuration) the ground-state energy 
is -28.30 MeV for Ml = 19 MeV and -27.56 MeV for 
Kl = 28 MeV, compared to the exact ^He CD-Bonn 
ground-state energy of —26.16 MeV The complete 
convergence of the ground-state energy can be obtained 
within the NCSM, as demonstrated, e.g., in Fig. 1 of 
Rcf. [l^. As expected, Fig. ^ shows that in the large 
model space the renormalization is weaker, i.e., there is 
less need for renormalization, so that the value obtained 
with the bare operator is similar to the value obtained 
with the renormalized operator. Thus, for 9^3 fm~^ 
one cannot distinguish between the results obtained us- 
ing bare operators in the 18/iri model space for HO fre- 
quency hn = 28 MeV and the ones obtained using ef- 
fective operators. For momenta q > 3 fm^^ the results 
using the bare operator deviate from the renormalized 
values, because the short-range correlations induced by 
the interactions are cast into the effective interaction. If 
one wants to account for short-range correlations using a 
bare short-range operator, one has to increase the model 
space, so that the effect of the short-range renormaliza- 
tion is negligible. However, such a scheme would require 
a vast number of hfl excitations to obtain a convergent 
result. 

In Fig. 121 we present pll{q) calculated in IQhil and 
IShft model spaces, with HO frequencies of 19 and 28 
MeV. In the upper panel we show the results obtained 
using bare operators. In this case, the values are spread 
over orders of magnitude. In contrast, the lower panel 
demonstrate independence of both model space and fre- 
quency, when using the appropriate effective operator. 
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FIG. 4: (Color online) Longitudinal-longitudinal distribution 
function for different model spaces and frequencies in ^^C, 
using bare (upper panel) and effective (lower panel) opera- 
tors. For the 4ftf2 model space, we show only the results 
using hil = 15 MeV. While the frequency dependence is not 
completely removed, the use of effective operators produces 
indistinguishable results at large q for different model spaces, 
as discussed in the text. 



although the ground-state energy is not converged, as il- 
lustrated in Fig. 13 Moreover, Fig. |21shows that the con- 
vergence depends strongly upon the HO frequency, when 
using bare operators. Thus, the results obtained with the 
bare operator in IShfl with hfl = 19 MeV are far from 
the results using effective operators; moreover, this curve 
shows a second minimum around q ~ 4.25 fm~^, whereas 
the converged results are almost flat and several orders of 
magnitude larger for this value of the momentum trans- 
fer. In contrast, even if still significantly different from 
the converged values, the results for hfl = 28 MeV are 
closer to the ones obtained with effective operators. 

One can better observe the influence of the frequency 
and model space in Fig. 0] where we present the results 
for ^^C. For ^^C, unlike the case of ^He, where we have 
used a Jacobi-coordinate HO basis (see, e.g., Ref. 0), 
the investigation was performed using the Many-Fermion 
Dynamics code |2C| , which employs a Slater determinant 
basis and becomes much more efficient for A > 5 than a 
Jacobi-coordinate approach. In this case, the number of 
many-body configurations increases very rapidly and one 
has to limit the truncation to a smaller model space. We 
present results for up to Ahfl model spaces. Again, note 
in the upper panel that the results obtained with the bare 
operators differ widely in shape and magnitude. When 
using effective operators, however, all curves collapse into 
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the same shape and agree with each other for g > 3 fm"-'^, 
as shown in the lower panel. Because the calculation is 
not fully converged, the minima still change significantly 
with frequency and model space, even when one uses ef- 
fective operators. One observes that, although the re- 
sults at high momentum transfer are very close together, 
a small dependence upon the HO frequency persists. 

In summary, we have investigated the longitudinal- 
longitudinal distribution function (part of the Coulomb 
sum rule) in the framework of the NCSM, utilizing the 
two-body cluster approximation. Thus, we have ex- 
tended our previou s ap plication of the effective operator 
formalism [ill Il3| | to the calculation of an observ- 
able that probes the short-range correlations. We find 
that even very small model spaces can provide an accu- 
rate description of the short-range observables, if effec- 
tive operators are employed. This investigation shows 
that reliable results can be obtained for short-range op- 
erators, even for heavier nuclei, such as ^^C, for which 



the Ofi.r2 results are accurate at higher q. As expected, 
intermediate- and long-range correlations can be best de- 
scribed by increasing the size of the model space and/or 
by a using higher order cluster approximation. 
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